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Abstract. We present the 2.4-45 ^ni ISO-SWS spec- 
trum of the Be star 7 Cas (BO. 5 IVe). The spectrum is 
characterised by a thermal continuum which can be well fit 
by a power-law Si, oc i/^-^^ over the entire SWS wavelength 
range. For an isothermal disc of ionized gas with con- 
stant opening angle, this correponds to a density gradient 



p(r) oc 



-2.8 



We report the detection of the Humphreys 



(6-00) bound-free jump in emission at 3.4 /zm. The size of 
the jump is sensitive to the electron temperature of the 
gas in the disc, and we find T f« 9 000 K, i.e. much lower 
than the stellar effective temperature (25 000-30000 K). 
The spectrum is dominated by numerous emission lines, 
mostly from Hi, but also some Hei lines are detected. 
Several spectral features cannot be identified. The line 
strengths of the Hi emission lines do not follow case B 
recombination line theory. The line strengths and widths 
suggest that many lines are optically thick and come from 
an inner, high density region with radius 3-5 R* and tem- 
perature above that of the bulk of the disc material. Only 
the a, /3 and 7 transitions of the series lines contain a 
contribution from the outer regions. The level populations 
deviate significantly from LTE and are highly influenced 
by the optically thick, local (disc) continuum radiation 
field. The inner disc may be rotating more rapidly than 
the stellar photosphere. 
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Be - Starsdndividual: 7 Cas - Stars: mass-loss - Infrared: 
stars 



Send offprint requests to: S. Hony (hony@astro.uva.nl) 

* based on observations obtained with ISO, an ESA project 
with i nstruments funded bv ESA Member states (especially 



1. Introduction 

Be stars are rapidly rotating main sequence or giant stars 
that are characterized by the presence of (variable) Ha 
emission, caused by high-density circumstellar gas. Often 
the Ha line profile is double-peaked and its width corre- 
lates with the projected rotational velocity (vsini) of the 
underlying star (e.g. Dachs et al. 1986 ). At infrared and 
radio wavelengths, the continuum energy distribution of 
Be stars is dominated by free- free and bound-free emission 
from the high-density circumstellar gas, which also causes 
the Ha emission. Direct imaging of the circumstellar ma- 



terial at radio wavelengths (Dougherty & Taylor 1995 ) and 
in the Ha line, (e.g. Stee et al. 199S) shows that the gas is 
in a disc-like geometry. This flattened geometry is also ev- 
ident from the optical continuum linear polarisation (e.g. 
Poeckcrt ct al. 197E) due to electron scattering. 



The physical mechanism responsible for the disc-like 
geometry is related to the rapid rotation of the star. Mod- 
els proposed in the literature include the wind compressed 
disc model (Bjorkman & CassincUi 1993), viscocity-driven 
outflow in the equatorial regions dOkazaki 19911 ) and non- 
radial pulsations ( Vogt et al. 1983 ). Which of these mod- 
els is correct can be tested by investigating the density 
and kinematical structure of the disc. We use the infrared 
spectral region to explore the structure of Be star discs. At 
these wavelengths disc emission dominates the spectrum, 
and a rich hydrogen and helium emission line spectrum 
is available. The infrared lines probe the inner regions of 
the disc, and their strength and width are important di- 
agnostic tools for the density and kinematical structure 
of the disc. Ground-based infrared spectra of 7 Cas were 



previously reported by e.g. Chalabaev & Maillard (1983), 
Lowe et a.1. (1985D and [Hamann fc Simon (19871). T he In- 
frared Space Observatory (ISO) ( pK^essler et al. 1996 ) with 
its Short W avelength Spectrometer (SWS) (De Graauw 



the PI countries: France, Germany, the Netherlands and the 
United Kingdom) with the participation of ISAS and NASA 



et al. 1996 ) is very well suited to explore the infrared part 



of the spectrum of Be stars. 
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In this study we present a preliminary analysis of the 
ISO-SWS spectrum of one of the brightest and best stud- 
ied Be stars in the sky, 7 Cas (B0.5IVe). We will show 
that the Hi emission line spectrum of 7 Cas is not well 
represented by Menzel case B recombination line theory. 
Many line fluxes correlate with the local continuum and 
are independent of the intrinsic line strength (Einstein A 
coefficient). The observed line fluxes and widths suggest 
that these lines are formed in an inner region with well- 
determined size, and that only the intrinsically strongest 
lines have a contribution from outer layers. This paper 
is organized as follows. In Sect. 2 we briefly discuss the 
observations and data reduction. Section 3 discusses the 
continuum and Sect. 4 deals with the line spectrum. Sec- 
tion 5 discusses some implications of our measurements 
for the structure of the disc of 7 Cas. 

2. Observations and data analysis 

The Be star 7 Cas was observed with the SWS on board 
ISO on July 22nd, 1996, as part of the SWS guaranteed 
time programme BESTARS. A full spectral scan (2.4-45 
/xm) using Astronomical Observation Template (AOT) no. 
1, speed 4 (De Graauw et al. 1996) was obtained, while 
also several AOT02 line scans were taken. The observa- 
tions were reduced using the SWS Interactive Analysis 
(lA^) software package, with calibration files equivalent to 
pipeline version 7.0. Further processing consisted of bad 
data removal and rebinning on an equidistant wavelength 
grid. The flux levels are accurate to within 5 per cent for 
the wavelengths shortward of 7 fim. The observations be- 
tween 7 and 12 /im (band 2C) suffer from memory effects; 
this has little influence on the measured line properties but 
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Fig. 2. Schematic representation of the method used to 
derive the temperature from the Humphreys jump. On 
the left the observed flux (drawn line) and the photo- 
spheric contribution (dotted line) are shown, on the right 
the normalized excess flux is shown. Also indicated is the 
wavelength (A') longward of the jump where the excess is 
equal to the excces at the jump. 



3. The continuum energy distribution 

The continuum energy distribution of 7 Cas at IR wave- 
lengths is dominated by free-free and bound-free emission 



from the ionized part of t h e circumstellar gas (e.g. Poeck- 



does ir crease the uncertainty of the continuum flux level ert fc Marlborough 197g ; Waters et al. 1987 ). The stel- 



to 15 per cent. At even longer wavelengths the signal to 
noise ratio decreases dramatically and only the strongest 
lines can be measured with reasonable accuracy. Most of 
the emission lines are partially resolved with a ratio of 
FWHM to the FWHM of the instrumental profile between 
1.4-3.5. Only 6 of the emission lines are considered unre- 
solved having this ratio below 1.4. Since the SWS instru- 



mental profile is approximately Gaussian (Valentijn et al 



lar contribution to the total flux is about 20 per cent at 
2.4 /im, based on extrapolation of a Kurucz model atmo- 
sphere fltted to the UV continuum ( Pelting et al. 1993 ). 
The spectrum can be well represented by a single power- 
law, Sjy oc I/", with a — 0.99 ± 0.05. This spectral slope 
is slightly, but significantly flatter, than that derived by 
Waters et al. (1987 ) , based on IRAS broad-band photome- 
try taken in 1983. We have used the simple isothermal disc 



1996) and the observed lines are well fitted by Gaussians, 



we estimate the original line width from: 



Wohs 



"^org 



2 , 2 



(1) 

is the original 

FWHM and Winst is the FWHM of the instrumental pro- 
e latter value varies with wavelength. To determine 



where 



is the observed FWHM, 



"^org 



ve use measured line widths of emission linos of 



file. Tb 

planetary nebulae; NGC 7027 and NGC 6302, observed in 
the same observing mode. No significant line profile vari- 
ations are observed. We show the final AOTOl spectrum 
in Fig. I. 



model of Waters (1986) to estimate the radial density gra- 
dient in the disc, assuming a power-law p(r) — po(r/R*)~", 
and find n — 2.8 ± 0.1. The value of po depends on the 
assumed opening angle 9 of the disc, as well as on the stel- 
lar radius and disc temperature. We use R* = 10 Rq and 
Tdisc = lO"' K (see below). Analysis of the optical linear 
polarisation and interferometric imaging of C Tau (Wood 



et al. 1997) suggests a half opening angle of 2.5°. We use 



a 1° half opening angle. The derived density at the stellar 
surface is po = 3.5x10"^^ g cm~^; an emission measure 
EM = 1.5x10*^^ cm~3 was found. 

There are some wavelength ranges that show a devi- 
ation from the power-law behavior of the continuum dis- 
cussed above. Near 3.28 /im the merging of the emission 
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Fig. 1. SWS AOTOl speed 4 spectrum of 7 Cas between 2 and 12 fim. The spectrum is dominated by numerous 
emission lines from hydrogen. A few He Unes are also observed. The inset shows the region aroud the Humphreys 
jump. The dashed line shows the jump in the continuum level due to the jump in bound— free opacity near 3.4 fim. 



lines of the Humphreys series, with lower quantum level 
n = 6, causes a jump. This Humphreys jump (seen in 
emission), which is similar to the Balmer jump at optical 
wavelengths, can be used to derive the average electron 
temperature of the emitting region. The difference in flux 
on both sides of the Humphreys jump is caused by a dis- 
continuity in bound-free (Kff+bf) opacity of the gas in the 
disc. We write for the total continuum opacity. 

Kg+M oc A" X (1 - e"=''/^'=^)/(c/i/Afcr) X r-3/2 x 

{g{X,T) + b{X,T)}, (2) 

where g(A,T) and b(A,T) are the free-free and bound-free 
gaunt factors, respectively. b(A,T) is a sensitive function 
of the temperature: the jump in b(A,T) (and in flux) in- 
creases towards lower electron temperature. The change 
in g(A ,T) is negligible over the wavelength range of in- 
terest. The jump in b(A,T) is thus a diagnostic of the 
temperature in the disc. We use the following method to 
determine the size of this jump: We define the normalized 
excess flux as Z^ — 1 = (Fa-Fa, *)/Fa^,, where Fa,* is the 
stellar photospheric flux, see Fig. ||. Za — 1 is normalized 
to the source function of the gas in the disc modulo a 



constant since both the disc and the star radiate in the 
Rayleigh-Jeans limit in this wavelength regime and thus 
have the same wavelength dependence. On the blue side 
of the discontinuity b(A,T) has a certain value, with a cor- 
responding value of Kff+bf , of Tff+bf for each line of sight 
through the disc and thus a corresponding value of Za — 1. 
Beyond the discontinuity there is a drop in b(A,T), Kg+bf , 
Tff-|_bf and Za — 1. Since there is a wide range of T^+bf for 
different lines of sight Za — 1 is not a simple function of 
b(A,T). However Kff+bf steadily increases with wavelength, 
e.g. Eqn. H, and thus there is a wavelength (A') where the 
loss in b(A,T) is compensated by the increase in A. At 
A' the Kff+bf is equal to the previous value, so Tff+bf and 
Za — 1 are also equal. We can determine A' directly from 
the observations, e.g. Fig. ||. 

Using A' = 3.470 ± 0.005 /im and the gaunt factors 



calculated by Waters & Lamers (1984) we find an electron 
temperature in the disc of 9 500 ± 1 000 K. This temper- 
ature would cause a weak but measurable jump at the 
Hansen-Strong series limit (near 4.5 /zm) but none is ob- 
served. However, we note that near this wavelength two 
strong H I lines are present that may mask an otherwise 
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detectable jump. The disc temperature agrees well with a 



density-weighted temperature of 10 700 K derived by Mil 
lar & Marlborough (1998) from an energy balance calcula- 
tion using the Poeckert & Marlborough model for the disc 
of 7 Cas. Note that the inner regions of the disc may have 
considerably higher temperatures, because this method is 
insensitive to contributions of those parts of the highest 
density parts of the disc where the continuum is optically 
thick, (see also Sect. 4). A second region which deviates 
from the power-law is near 4.3 fj,m. We cannot find a rea- 
sonable identification for this spectral feature. The CO2 
stretch mode band is at 4.27 fj,m. However, the interstellar 
extinction towards 7 Cas is very low which implies that 
we must rule out this possibility. 



4. The emission lines 

The entire SWS spectral region, but especially bands 1 
and 2 (2.4-12 /im), is dominated by strong and partially 
resolved emission lines. The vast majority of these lines 
are Hi recombination lines. We find the series limit of 
the Humphreys and Hansen-Strong series, as well as those 
from lower levels 8, 9 and 10. Several lines with lower lev- 
els above 10 were also identified. We find a few Hei lines 
(notably at 2.4861, 2.5729, 4.2960 and 4.0367 nm.). No 
forbidden lines could be found, although several emission 
features are unidentified and could perhaps be attributed 
to forbidden lines. The strongest unidentified line is at 
2.8934 /zm, close to the 2.8964 [Niii] fine. The lack of fine- 
structure lines is consistent with the optical spectrum of 
7 Cas. It is likely that the high density of the ionized gas 
causes collisional de-excitation of the fine-structure tran- 
sitions. This situation is markedly different for the hy- 
pergiant P Cygni, whose infrared spectrum is also domi- 
nated by Hi recombination lines from circumstellar gas, 
but which also shows prominent emission fr om e.g. [Fe 11 



[Niii] and [Si 11] in its ISO-SWS spectrum ( [Lamers et al 



1996) 



We have measured the line fluxes with respect to both 
the local continuum and the stellar photospheric contin- 
uum at line centre. The measured line properties are given 
in Table |l|. In Fig. ^ we show the resulting curve of growth 
for the emission lines, where we plot the equivalent width 
EW divided by wavelength versus log Xuno- This last 
quantity is proportional to the line optical depth ( ^aal 
et al. 1995| ). The observed line fluxes (Fig. ||) are all much 
smaller (a factor 10 or more) than expected on the ba- 
sis of the continuum emission measure determination and 
the optically thin case B recombination line predictions 



(Hummer & Storey 1987). This shows that even for the 



weakest lines the emission is optically thick: the combined 
line and continuum opacity is much larger than 1. This 
is not surpri sing given the shape of the continuum energy 
distribution; Waters et al. (1991 ) show that the continuum 
becomes optically thick near 0.8 /im. 
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Fig. 3. Comparison of measured line strengths with pre- 
dictions based on optically thin case B recombination the- 
ory. The dotted line denotes the locus where the measured 
fluxes would equal the predictions. The crosses correspond 
to Brackett a and /3, the open squares to the Pfund se- 
ries, the filled squares to the Humphreys series, the open 
triangles to the Hansen-Strong series, the stars to n = 8, 
and the filled triangles to n = 9. 



The shape of the empirical curve of growth is very 
different from that expected on the basis of LTE line for- 
mation in a rotating, partially optically thick disc. Model 
calculations show that under these conditions the curve of 
growth has a linear part (where line emission is optically 
thin and proportional to line strength), and a power-law 
part whose slope depends on the radial density gradient of 
the gas ( |Zaal et al. 1995[ ), see also [Wellmann (1952| ). The 
observed curve of growth however shows a very steep rise 
of EW/A with Xiinc up to log Xunc ~ —36.5, followed by 
a flat part for —36.5 < log Xuno < —35.0 whose slope is 
close to zero, i.e. much smaller than 0.4 expected on the 
basis of the density gradient derived from the continuum 
free-free and bound-free excess. Finally, a rising part for 
the strongest lines in the spectrum (log Xunc > —35.0) is 
seen. The flat part of the curve of growth is surprising and 
suggests that these lines have saturated and are formed in 
a region with a well-defined outer radius. We note that 
the size of this region is not the same for every series, but 
increases with lower quantum level. 

In order to understand better the nature of the line 
formation in 7 Cas, we show in Fig. |^ the line strength 
EW/X, where EW is the line equivalent width measured 
with respect to the local continuum, as a function of wave- 
length. Also shown in Fig. || is the line FWHM versus 
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Table 1. Properties of hydrogenic emission lines. A is the peak wavelength of the Gaussian fit. Typical errors on these 
are l/2500th of the wavelength, width is the FWHM of the line after deconvolution. / is the integrated line-flux. cont. 
is the mean continuum level underneath the line. The errors on the continuum level are dominated by the absolute 
flux calibration uncertainty of the SWS instrument. 
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wavelength. Both EW /\ and the FWHM show a char- 
acteristic wavelength dependence: for each series, EW /\ 
increases and then saturates to a value of 5-6x10"^, while 
the line FWHM decreases and reaches a roughly constant 
value of 250-300 km s~^. In other words, the lines reach a 
constant line over continuum ratio of about 1.3 at a con- 



stant FWHM of 250-300 km s'^ Only the q: and /3 lines 
of each series deviate from this behaviour: their line over 
continuum ratio is larger and their FWHM smaller. 

The decrease of FWHM with increasing line strength 

is expected for lines formed in a rotating disc in which 
the rotational velocity decreases with distance; as more 
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Fig. 5. Top panel: EW/A of the Hi lines in the ISO-SWS spectrum of 7 Cas as a function of wavelength. Symbols 
are the same as in Fig. ||. The EWs have been determined with respect to the local continuum. The EWs for all series 
behave in a similar way. For moderately strong lines the EW/A does not depend on wavelength and is about 5x10"'^ 
irrespective of the series. Only the strongest lines of each series deviate from this trend. Bottom panel: observed full 
width at half maximum of the H i lines versus wavelength. The symbols refer to the same series as in the upper panel. 
The lines become narrower with increasing line strength as expected in a rapidly rotating disc in which the rotational 
velocity decreases with distance. The FWHM does not decrease further for moderately strong lines. These lines also 
show a constant EW/A (upper panel). 



line flux is coming from the outer, more slowly rotating 
regions, the line width will decrease. This is direct proof 
of the rotating nature of the line emitting region. It is 
remarkable that the weakest lines of each series have a 
very high FWHM of more than 550 km s~^. Such high 
velocities are not expected given the photospheric vsini 



of 230 km s-i (Slettebak 1982). This suggests that the 



inner disc is rotating more rapidly than the star. However, 
line broadening due to electron scattering could also cause 
these large line widths, but we do not observe prominent 
electron scattering wings. 

If we assume that the line and continuum source func- 
tions are equal, no line emission from the layers with 
Tfi > 1 should be detectable. However, the large line width 
of the weakest lines strongly suggests that this line emis- 
sion is originating from rapidly rotating parts of the disc. 
The fact that the line FWHM decreases with increasing 



line strength shows that the broad, weak lines are formed 
closest to the star. In the 2.5 to 7 ^m wavelength region, 
where most of the weakest lines in Fig. |^ are located, the 
continuum is optically thick out to a radius of 2.5 to 3 R*. 
Assuming Keplerian rotation in the disc, the part of the 
disc which is not optically thick for continuum radiation 
rotates at projected speeds less than about vosini/1.7, 
where vq is the Keplerian velocity of the disc at the stel- 
lar surface. Using reasonable values for the mass and ra- 
dius for a B0.5 IV star of R^^IORq and M*=15Mq , 
we find a maximum speed of about 220 km s^^. Emis- 
sion lines in the disc whose source function is equal to 
that of the continuum therefore should have FWZI less 
than 440 km s^^, but even the FWHM of the weak fines 
is significantly higher than this. It is unlikely that rota- 
tional velocities are as high as 250-300 km s~^ as far out 
as 2.5-3 R*, unless the rotational velocity field deviates 
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Fig. 4. Measured line fluxes, expressed in equivalent width 
(EW/A) with respect to the stellar continuum versus line 
opacity. The line fluxes are independent of line strength 
between log Xuno ~ —36 and —35. 



significantly from Keplerian. We conclude that we detect 
line emission from those parts of the disc that are opti- 
cally thick for continuum radiation, and hence the line 
source function in these inner regions must exceed that 
of the continuum (which is the Planck function at the lo- 
cal electron temperature). Several effects can cause such 
a larger source function: NLTE level populations, or an 
enhanced electron temperature near the upper part of the 
disc. If the latter effect is important, temperatures in the 
line forming regions may be as much as 30 per cent higher 
than in the bulk of the disc, since the stronger lines are 
1.3 times the continuum. 

As pointed out above, the observed line fluxes of every 
series first increase with intrinsic line strength, but then 
reach a constant value with respect to the local contin- 
uum, both in EW/X and in the line to continuum ratio. 
This effect occurs for transitions from upper levels roughly 
between 15 and about 10, to lower levels between 9 and 5. 
The line fluxes are no longer determined by the intrinsic 
line strength, but by the local continuum! The contribu- 
tion from outer regions, with lower rotational velocities, 
is absent. This can be seen from the FWHM's of the lines 
with constant line over continuum ratio, that no longer 
decrease but are constant at about 250-300 km s~^. This 
could be due to a finite outer radius of the disc. However, 
the ISO-SWS spectrum as well as ISO-PHOT photometry 
at 90 /Ltm (Trams et al., in preparation) show that there 
is no change in continuum slope over a wide wavelength 
range. Such a change in slope would be expected in the 



case of a finite outer radius. We conclude that the lack of 
line emission from the outer regions cannot be due to a 
finite outer disc radius. We note that the fact that the line 
to continuum ratio is constant for certain lines, irrespec- 
tive of the strength of the local continuum and thus of the 
size of the continuum emitting region, also argues against 
a finite outer radius. 

The lack of line flux from the outer regions therefore 
must be due to a change (decrease) in the line source func- 
tion compared to the inner regions. For a certain range of 
intermediate energy levels in the atom, the optically thick 
IR free-free disc radiation field in the inner regions is able 
to maintain significant level populations, but these lev- 
els rapidly become de-populated when the local free-free 
continuum becomes optically thin. Therefore the line flux 
does not increase with intrinsic line strength (as do the 
weaker lines from high upper levels), but reaches a line 
over continuum ratio determined by the ratio of source 
functions of line and continuum. 

5. Discussion 



Hamann & Simon (1987) discuss the near-IR spectrum of 
7 Gas and note that the line flux ratios of the weak Pfund 
lines to the Br7 line are very difficult to reconcile with 
those of Br7 and Bra for any reasonable choice of radial 
density gradient in the disc. We confirm this result and 
show that it is a general property of the line fluxes of 7 Gas 
to behave in a highly non-standard way, probably caused 
by a combination of non-LTE and temperature effects. 
We do not confirm the lack of correlation between line 



strength and line width noted by Hamann fc Simon (1987 ); 
the ISO data clearly show such a correlation for all series 
in our spectrum. This may be due to the Earth atmosphere 
affecting the estimate of the weak high level Pfund series 
lines in the spectrum of Hamann & Simon. The line widths 
of their stronger lines (e.g. Bra) agree well with those of 
our spectrum. They also agree well with the line widths 
derived from high-resolution spectra published by Lowe 
|et al. (1985| ) and phalabaev fc Maillard (1983] ). 

The FWHM of the weakest H i lines (more than 500 
km s^^) significantly exceeds that of the photospheric 
vsini of 230 km s"^. While electron scattering may af- 
fect the width of the strongest lines (especially at short 
wavelengths where free-free opacities are small compared 
to the electron scattering opacity) the IR lines do not 
show evidence for prominent electron scattering wings. 
This suggests that the line width is mainly due to kine- 
matical broadening, which would imply that the disc is 
rotating more rapidly than the underlying star (ignoring 
errors in v sin i) . This can only occur if some mechanism 
adds angular momentum to the material injected into the 



disc. We note that Hanuschik (1987) found line widths 



of weak Fell emission lines in the optical spectra of Be 
stars that also significantly exceeded the width expected 
on the basis of the observed vsini. Transfer of angular mo- 
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mentum may lead to spin-down of the star ( Porter 1998| ). 
Using the formahsm given by Porter (199^ ), and assum- 
ing a full opening angle of 2 degrees and po of 3xl0~^^ 
g cm~'^, the outflow velocity in the disc near the star can 
be of the order of 1 km without significant spin-down 
of the star during its main sequence life time. Such a value 
is in agreement with the observed line shape. 

The picture which emerges from our analysis of the in- 
frared spectrum of 7 Cas is that of a circumstellar region 
of very high density, perhaps exceeding 3 10~^^ g cm~^, 
which is rotating rapidly and whose rotational velocity de- 
creases with distance. The rotational velocity in the disc 
near the star exceeds that of the photosphere. This re- 
gion is heated by radiation from the central star, and the 
surface layers which directly absorb the stellar radiation 
field have higher temperatures than regions closer to the 
equatorial plane of the disc. While the weak lines originate 
from these dense, warm regions, only the strongest, a and 
/?, lines of each series probe the outer portions of the disc. 
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